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A beam's density matrix that is described by the superposition of excitation on coherent states 
with thermal noise (SECST) is presented, and its matrix elements in Fock space are calculated. The 
maximum information transmitted by the SECST beam is derived. It is more than that by coherent 
light beam and increases as the excitation photon number increases. In addition, the nonclassicality 
of density matrix is demonstrated by calculating its Wigner function. 
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I. INTRODUCTION 

Recently, much attention has been paid to the exci- 
tation on coherent states (ECS) P, [3, i, i, i, i]. As 
pointed out in Rcfs.@,@|) the single photon ECS causes a 
classical-to-quantum (nonclassical) transition. The ECSs 
can be considered as a generalization of coherent states 
P, Q and number eigenstates. All these states can be 
used as signal beams in optical communications field, in 
which the nonclassicality of signals plays an important 
role. 

However, in reality, signal beams are usually mixed 
with thermal noise. Statistical properties of the super- 
position of (squeezed) coherent states with thermal light 
(SCST) have been investigated by calculating the pho- 
ton number matrix elements (A^| p \M) of SCST's density 
matrix [9|, llO| . These properties are useful in quantum 
optics and quantum electronics (e.g. how lasers work- 
ing well above threshold, heterodyne detection of light, 
etc.) [11| . Some general properties of the density matri- 
ces which describe coherent, squeezed and number eigen- 
states in thermal noise are studied in Ref. [l2| • It is found 
that the information transmitted by the superposition of 
number eigenstates with thermal light (SNET) beam is 
less than that by the SCST beam [13|. 

In this paper, we investigate statistical properties of 
the superposition of ECS with thermal hght (SECST). 
We present the relevant density matrix in Fock space 
and derive the Mandel Q parameter. The SECST field 
can exhibit a significant amount of super-Poissonian pho- 
ton statistics (PPS) due to the presence of thermal noise 
for excitation photon number m — Q\ while for to ^ 
the SECST field can present the sub-PPS when the ther- 
mal mean photon number is less than a threshold value. 
In addition, the threshold value increases as m increases. 
We also calculate the maximum information (channel ca- 
pacity) transmitted by the SECST beam, which increases 
as TO. increases. In addition, the nonclassicality of den- 
sity matrix is also presented by calculating the Wigner 
function of the SECST. 

Our paper is arranged as follows. In Sec. II we present 
the density matrix p that describes the SECST and calcu- 
late its matrix elements in Fock space by using the normal 



ordered form of p. The PPS distributions are discussed in 
Sec III. The maximum information is calculated in Sec. 
IV. Sec. V is devoted to deriving the Wigner function of 
the SECST and discussing its nonclassicality in details. 
Conclusions are summarized in the last section. 



II. EXCITATION ON COHERENT STATES 
WITH THERMAL NOISE 

Firstly, let us briefly review the excitation on coherent 
states (ECSs). The ECSs, first introduced by Agarwal 
and Tara [if , are the result of successive elementary one- 
photon excitations of a coherent state, and is an inter- 
mediate state in between the Fock state and the coher- 
ent state, since it exhibits the sub-Poissonian character. 
Theoretically, the ECSs can be obtained by repeatedly 
operating the photon creation operator a^ on a coherent 
state, so its density operator is 



Po = a,ma^™ |«) («| a" 



(1) 



where Ca,m = ["^!^m(~ |q^| )] ^ is the normalization fac- 
tor |a) = exp(— |a| /2 -I- aa^) |0) is the coherent state 
[71, |8| , and Lm {x) is the mth-order Laguerre polynomial. 
The SECST is described by the density matrix [12] 
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where D (z) — exp(zat — z*a) is the displacement oper- 
ator, and ht is the mean number of thermal photons for 
Pa -^ |0) (0|- We can easily prove that Trp — 1, as it 
should be. In fact. 
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A. Normal ordering form of the SECST 

For the simplicity in our later calculation, we first per- 
form the integration in Eq.([2|) by using the technique of 
integration within an ordered product (IWOP) of oper- 
ators [IJj \lM- Using the normal ordering form of the 
vacuum projector |0) (0| =: exp(— a^a): , we can reform 
Eq.([2]) as the following form 
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(5) 



X exp [(a^ — a*) z + {a — a) z*] 
X [a' ^ z j [a — z) : . 

In the last step of (O, we noticed that for any operator 
/(at, a) 

D [z) f{a\a)D^ (z) = f{a^ ~z*,a~ z). (6) 

Making two independent variable displacements, 



r,a- 



/3, 



(note that operators a^ , a can be considered as C-number 
within the normal order : : ), thus Eq.([5]) can be rewrit- 
ten as 
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where we have set \t = ^Jnt/{^- + nt) and A = \t{^a + 
a). Then using the integration expression of two- variable 
Hermite polynomial Hm,n [16l |. 
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we can put Eq.([7]) into 
P 



C^^2m+2._ i^iY^H,^.,^{AK-A) 



nt 



X exp 
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nt + 1 



(9) 



In particular, when m = 0, corresponding to the case of 
superposition of coherent state with thermal noise, Eq.Q 
reduces to 



1 



nt + l 



D{a) : e ^TTT ■, ^t (q,) ^ 



(10) 



which can be directly checked by using Eqs.([2) and ^ 
as well as noticing po = \a) {a\ . 

Further employing the relation between Hermite poly- 
nomial and Laguerre polynomial 1(^, 



tl'm,n (,Sj ^) 



we can see that 
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(12) 



Eqs.© and (|12p are the normal ordering form of the 
SECST. From these it is convenient to calculate the phase 
space distributions, such as Q-function, P-representation 
and Wigner function. 



B. The matrix elements {N\ p\M) 

Now we calculate the matrix elements of p in Eq.Q 
between two number states {N\ and \M) , i.e., (iV| p \M) . 
Employing the overcompleteness of coherent states, one 
can express the matrix elements (A^| p \M) as 

(iV| p \M) ^ I ^^ {N 1/3) (/3| p I7) (7 |M) , (13) 

where the overlap between the coherent state and the 
number state is given by 

1 



(7 \M) 



'M\ 



e-l7lV2^*A/^ 



(14) 



and the matrix elements (/3| p I7) can be obtained from 
Eq.(l9]) due to p's normal ordering form. 
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where we have used the generating function of two- 
variable Hermite polynomial Hm.m 



H,n,n (X, y) = g^,rig^fn °^P t"^*' + ^^ + ^'^1 



(16) 



When M = N, {N\ p \N) is just the photon number dis- 
tribution of the SECST. Then combing with Egs. pS)) . 
(fT3| and (fT4| . after a lengthy but straightforward calcu- 
lation, one can get the matrix elements {N\ p \M) , (with- 
out loss of generality, let M ^ N) 
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where we have used the integral formula |17J 
^/(/3*)exp{-|/3|Vr/3}=/(T), 



(17) 



(18) 



and another expression of two-variable Hermite polyno- 
mial Hm,n, 

min(m,n) , ^ .1 ^m-l n-l 
H^.n{£.,K)= 2_^ —^^j- j^j-j— j^^. (19) 
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In particular, when m = 0, noticing M ^ N and Eq. pTj) . 
Eq.([T7l) reduces to 
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, (20) 



which is just the Glauber-Lachs formula [9| when fit — 
(e^" — 1)~^. While for a = 0, corresponding to the case of 
superposition of number state with thermal light, using 
Eq.dlSl), Eq.dni) becomes 



{N\p\M) ^Sm,nPn, 
where {kg — max[0, m — N]) 

k+N 



P, 



iim 
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nt 
1 lfit + 1 



[fit {fit + 1)] 



k—r 



nt + l ^ fc! [k + N ~ my. [(m - 
Eq.dH]) is just the result of Ref. [13]. 
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III. SUB-POISSONIAN PHOTON STATISTICS 

To see clearly the photon statistics properties of the 
SECST, in this section, we pay our attention to the 

variance of the photon number operator ( (An) ) — 

(n^) — (n) ■ In particular, we will examine the evolu- 
tion of the Mandcl Q parameter defined as 

Q = A^-^ - (at a) 
{a'a} 

_ (a^at^) - (aat)' - (aa^) 



it)-l 



(23) 



which measures the derivation of the variance of the pho- 
ton number distribution of the field state under consid- 
eration from the Poissonian distribution of the coherent 
state. Q = 1,Q > 1 and Q < 1 correspond to Poisso- 
nian photon statistics (PPS), supcr-PPS and sub-PPS, 
respectively. 

In order to calculate the average value in Eq. ((23|) . we 
first calculate the value of {a\ a^a''"^ \a) . In fact, using 



ra+n \m+n i \ 



a 



\a) = {m + ny.Lm+n{-H) (24) 



and 



I z"z*"'P(z) = n1'm\6r, 



(25) 



we can evaluate (for writing's convenience, let L„i denote 

Lm{~\a\^)) 
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(26) 



and 



(a^a^^) ^ 2nl + — [mtL^+i + (m -f 2) L^+2] 

Substituting Eqs.dSB]) and ^ into ^ leads to 
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(28) 



At the zero-temperature limit {nt — > 0), Eq. ((28|) just 
reduces to Eq.(2.20) in Ref.[l|. 

In Fig.l, we display the parameter Q {nt, \a\) as a func- 
tion of {nt, I a I) for different values m. From Fig.l, we see 
that, for the excitation photon number m = (see Fig.l 
(a)), Q {fit = 0, \a\) = 1 corresponding to coherent state 
(a PPS); while Q {fit ^ 0, |a|) > 1, i.e., the SECST field 
exhibits a significant amount of super-PPS due to the 
presence of n*. From Fig.l (b) and (c), we see that, when 
m ^ 0, the SECST field presents the sub-PPS when fit 
is less than a threshold value for a given |a| ; the thresh- 
old value increases as m increases. For example, when 
\a\ = 0, the threshold values are about 0.414 and 0.481, 
respectively, for m — 1 and m — 6. 



IV. INFORMATION TRANSMITTED BY THE 
SECST BEAM 

According to the negentropy principle of Brillouin [i^| , 
the maximum information / transmitted by a beam is 



I^Sn 



Jact, 



(29) 



in which S'max and Sact represent the maximum entropy 
and the actual entropy, respectively, possessed by the 



QlSif.ti-) 





FIG. 1: (Color online) The evolution of Mandel Q parameter 
as a function of (wt, |a|) for different values m. 



quantum mechanical system described by a density ma- 
trix p. Here the maximum information / is an ideal one 
transmitted through an ideal optical communication sys- 
tem. 

For the SECST system, the actual entropy is 



bac 



-Tr(plnp) = -y^o-Arlno-AT, 



(30) 



N 



where p = J2n '^n \N) {N\ , and ctat — {N\ p \N) . (Jn can 
be obtained from Eq. fTT)) . i.e., 
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w— a.t; — a 



which is independent of the phase of a. On the other 
hand, for a system in thermal equilibrium, described by 
the density matrix ptin with mean photons number nt. 





FIG. 2: (Color online) The maximum information I (nt, |Q|)as 
a function of {fit, \a\) for some different values (a) m = Q, (6) 
m — 1, (truncating the infinite sum at A^max ~ 70). 



its entropy is 



S" = - ^ Pat In Ptv = hr(l + fit) + fit In 



N 



fit + l 
nt 



, (32) 



sN=l 



where Pjv = n^ / {nt + 1)"~^ obtained from Eg. ipOl) un- 
der the condition m — 0,a — 0. Note that the maximum 
entropy of the system is equal to the entropy of a system 
in thermal equilibrium, with an equal mean number of 
photons. The mean photons number of the SECST is 
given by Eq. ipS)) . Therefore, using Eg. ([5^ . we have 
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nt 
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(1 + to) Lm+i + ntLm 
(H- m) L„i+i + (fit - 1) Pr; 



(33) 



From Eqs. ([29|) . ([30|) and ([33|) . we can calculate the 
maximum information transmitted by the SECST beam. 
In Fig. 2, the maximum information I {fit, \a\) is plot- 
ted as a function of {fit, \a\) for some different values m 
(truncating the infinite sum at iVmax — 70). From Fig. 2, 
we can see that, for a given fit, I {fit, |a|) increases as the 
value |a| increases; for given |a| , in general, 7 (fit, |q|) 
grows up as the value fit increases. In order to see 
clearly the effect of different parameter m to I {fit, |a|) , 
we presented a plot in Fig. 3, from which it is obvious 




FIG. 3: (Color online) The maximum information 
I {nt,\a\ = l)as a function of (fit) for some different values 
(a) m — 0,{b) m = 1, (c) m = 2 (truncating the infinite sum 
at iVmax = 70). 



that / (fit, \a\) becomes bigger due to the presence of m, 
and increases as m increases. In other words, the max- 
imum information transmitted by the SECST beam is 
larger than that by the SCST (?ti = 0). The channel of 
ECS can carry with more information than that of co- 
herent state. In Ref. [ij], Vourdas has pointed out that 
the coherent signals (of known phase) can transmit more 
information than the number eigenvectors signals. Thus 
among these three beams, the SECST beam can transmit 
most information. 
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THE WIGNER FUNCTION OF THE SECST 



A. The Wigner function 



The Wigner function (WF) plays an important role in 
quantum optics, especially the WF can be reconstructed 
from measurements [1^ [20| • The WF is a powerful tool 
to investigate the nonclassicality of optical fields. The 
presence of negativity in the WF of optical field is a sig- 
nature of its nonclassicality is often used to describe the 
decoherence of quantum states. In this section, using 
the normally ordered form of the SECST, we evaluate its 
WF. For a single-mode system, the WF is given by [21| 



1^(7,7*) 



=2|7|- 



d2/3 



-/3|p|/3) 



„2(/3*7-/37*) 



(34) 



where |/3) is the coherent state and 7 = x + iy. From 
Ea. ([M)l it is easy to see that once the normal ordered 
form of p is known, we can conveniently obtain the WF 
of p. 

On substituting Eq.^ into Eq.(I331) we obtain the WF 



FIG. 4: (Color online) The evolution of the Wigner function 
of the SECST with a = 0.2 -I- 0.2i for several different values 
m and fit. 



of the SECST 
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where we have set 
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Noticing that L„^{- \a\^) > 0, and i„[- |A2|^/Af] > 



when 1 



> 0, thus the WF of the SECST 



(2rH + l)nt 

is always positive under the condition of fif > 1/2. In 



particular, when m — 0, Eq. (j35p becomes 
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(37) 



which corresponds to the 
mean photon number fit. 

4|7|'(n, + l)/{(2n, + l)[(2nH 
([55| yields 



thermal state with 
While for a = 0, 
l)n,], \A2t/Al ^ 
l)nt-l\} = e, Eq. 
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[{2nt + l)nt-lX 



^im(-0 
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(38) 
At the zero-temperature limit, T ^ 0,nf ^ 0, Eq. ((37|) 
reduces now into - exp(— 2 \a — 7I ), i.e., the WF of co- 
herent state (a Guassian form), which can be seen from 
Eq.© yielding p = \a) {a\ under the condition m = Q; 
while Eq.® becomes i(-l)™e-2|^l'L„(4 |7|^), corre- 
sponding to the WF of number state. 

Using Eq.dini), the WFs of the SECST are depicted 
in Fig. 4 in phase space with a = 0.2 -I- 0.2z for several 
different values m and fit . It is easy to see that the nega- 
tive region of WF gradually disappears as m and fit and 
increases. 



B. The marginal distributions of the SECST 

We now find the probability distribution of position or 

momentum! the marginal distributions, by performing 

the WF either over the variable y or the variable x, re- 
spectively. Using Eqs. ((35)) and ((36)) we can derive (denote 

J = X + iy, a = q + ip) 



P{x,nt) = / W{x,y)dy 



v/27r(2nt + l) (2^* - 1)' 
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■)2fe-r 



\H, 



rn — k 



(i^i)r 



(39) 



where Hm {x) is single variable Hermite polynomial and 
El = [{2fit - 1) a -t- 2x -I- 2ip] / ^2 (2nt + 1). Eq.® is 
the marginal distribution of WF of the SECST in "x- 
direction" . 



On the other hand, performing the integration over dx 
yields the other marginal distribution in "y-direction" , 



P (2/, nt) 



i^uirc^. 



V'27r (2nt -h 1) (2nt - 1)' 






\Hm-k{E2)\\ (40) 



where £'2 = i {2fita -a + 2q + 2iy) /[^J2{2fit + 1)]. As 
expected, the two marginal distributions are both real. 

VI. CONCLUSIONS 

In summary, we have investigated the photon statistics 
properties of the SECST, described by the density matrix 
p ^. We have calculated the matrix elements (A^| p \M) 
in Fock space and the Mandel Q parameter. It is found 
that the SECST field exhibits a significant amount of 
super-PPS due to the presence of thermal noise (fit) for 
excitation photon number m — Q and that, for m ^ Q and 
a given |a| , the SECST field presents the sub-PPS when 
fit is less than a threshold value. In addition, the thresh- 
old value increases as m increases. We have presented 
the maximum information (channel capacity) transmit- 
ted by the SECST beam. It is shown that the maxi- 
mum information transmitted increases as m increases. 
This implies that among the coherent signals, the eigen- 
number signals and the ECS in thermal light, the last 
one can transmit the most information. Further, as one 
of the photon statistical properties, the Wigner function 
and the marginal distributions of the SECST have also 
been derived, from which one can clearly see the nonclas- 
sicality. The negative region has no chance to be present 
when the average photon number fit of thermal noise ex- 
ceeds 1/2. The marginal distributions are related to the 
Hermite polynomial. 
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